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Local correlation functional for electrons in two dimensions 
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We derive a local approximation for the correlation energy in two-dimensional electronic systems. 
In the derivation we follow the scheme originally developed by CoUe and Salvetti for three dimen- 
sions, and consider a Gaussian approximation for the pair density. Then, we introduce an ad-hoc 
modification which better accounts for both the long-range correlation, and the kinetic-energy con- 
tribution to the correlation energy. The resulting functional is local, and depends parametrically on 
the number of electrons in the system. We apply this functional to the homogeneous electron gas and 
to a set of two-dimensional quantum dots covering a wide range of electron densities and thus various 
amounts of correlation. In all test cases we find an excellent agreement between our results and the 
exact correlation energies. Our correlation functional has a form that is simple and straightforward 
to implement, but broadly outperforms the commonly used local-density approximation. 
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I. INTRODUCTION 

In the last couple of decades, the growing world of nan- 
otechnology put at our disposal several classes of low- 
dimensional materials. Particularly interesting exam- 
ples are two-dimensional (2D) quantum dotsA'- (QDs), 
formed at the interface between two semiconductors. 
These systems are not only important from a technolog- 
ical point of view, but are also remarkable from a purely 
theoretical perspective. In fact, as they can be built with 
different shapes and sizes, and with a varying number of 
electrons, they are the ideal system to study electronic 
correlation. 

The problem of electronic correlation is perhaps the 
most challenging in the field of condensed matter physics. 
Numerous approaches to handle this problem, with vary- 
ing degrees of sophistication and complexity, have been 
put forward since the very birth of quantum mechan- 
ics. Few-electron QDs can be studied accurately by, e.g., 
configuration interaction^ (CI), or by quantum Monte 
Carlo techniques^'^'^ (QMC). To describe the electronic 
properties of larger dots one has to resort to alterna- 
tive approaches such as extended Hartree-Fock^ (HF) or 
density- functional theory^i^i^ (DFT). 

In DFT, the complexity of the many-body problem 
is embodied in the so-called exchange and correlation 
functional. Several approximations exist for this quan- 
tity, allowing for very accurate calculations of electronic 
properties in atoms, molecules, and solids. Clearly, most 
exchange-correlation functionals are derived for three- 
dimensional electronic systems. However, these approx- 
imations are known to break down when applied in the 
2D limit ji^ This calls for new formulas specialized for the 
2D case. Particularly challenging in these applications is 
the fact that, compared with atomic systems, correlation 
effects in 2D typically have a more prominent role due to 
the large size of the systems (from 10~* to 10~^ m), and 



to their low electronic densities. 

Within DFT, 2D systems such as QDs are commonly 
studied using the 2D version of the local-density approxi- 
mation (LDA) . It is a combination of the exchange func- 
tional derived for the uniform 2D electron gas by Ra- 
jagopal and Kimball^^- and the corresponding correla- 
tion functional fitted to accurate QMC calculations. The 
first of these LDA correlation functionals was put for- 
ward by Tanatar and Ceperleyi^ in 1989. Later on, it 
was generalized for the complete range of collinear spin 
polarizations by Attaccalite et alr^ Applications of the 
2D-LDA to QDs have been generally successful, even up 
to high magnetic fieldsi^i^ The LDA, however, suffers 
from several shortcomings, already well known from the 
three-dimensional world, especially for strongly inhomo- 
geneous systems, or in the low-density (strong correla- 
tion) regime. 

Several alternative paths exist to go beyond the simple 
LDA. A particularly successful approach starts with the 
seminal work of CoUe and Salvettii^ii^ (CS) who, start- 
ing with a physically motivated ansatz for the many-body 
wavefunction, developed a closed formula for the corre- 
lation energy. This formula has received a large inter- 
est, especially because it was used to derive the popular 
Lee- Yang-Parr (LYP) generalized gradient functional: — 
Together with Becke's exchange functional^'* it forms the 
BLYP functional, and in hybrid schemes it is a part of 
B3LYP^ X3LYP^ etc. 

Interestingly, the same CS formula can also be inter- 
preted as an orbital-dependent correlation functional, es- 
pecially suited for DFT calculations beyond the exact- 
exchange.-^^ It should, however, be emphasized that the 
CS correlation-energy functional has several known lim- 
itationsi^^iS^i^ In particular, while short-range correla- 
tions are well described,— important long-range correla- 
tions are missing. Even if these latter effects often can- 
not be ignored in large molecules and solids, they can be 



energetically negligible in small systems such as atoms. 
However, it has been shown recently that the long-range 
correlation problem may be cured to some extent i^^i^^ 
Secondly, in the CS functional the kinetic-energy con- 
tribution to the correlation energy (named below as the 
kinetic-energy correlation) is taken into account only in 
an empirical fashion through the fitting parameter. In 
this context, an interesting modifications of the original 
CS approach have been recently proposedi^ 

In this work, we generalize the CS scheme^»i^ to 2D. 
Then we use a Gaussian approximation for the pair prob- 
ability function. Finally, we introduce an ad-hoc modi- 
fication which, post-factum, seems to recover both the 
long-range and the kinetic-energy correlation to some 
good extent. 



II. THEORY 

Our starting point is the following ansat a^^'^^ for the 
many-body wavefunction ^P 

*(ricri,...,r7vcr7v) = *SD(riCri,...,r AT crjv) 

x[][l-^(r„r,)] . (1) 

Here, r and a denote respectively the space and spin 
coordinates of the electrons, and ^sd indicates the sin- 
gle Slater determinant of HF theory, which in the DFT 
context should be replaced by the Slater determinant 
generated from the occupied Kohn-Sham orbitals. The 
function (p describes the correlated part of the wavefunc- 
tion. In the center-of-mass, r = (ri -|-r2)/2, and relative, 
s = ri — r2, coordinate system, it can be written as 



(p(r,s) = [1 -$(r)(l + as) 



-I3\r)s' 



(2) 



where the quantities <&, a, and /3 act as correlation fac- 
tors. We point out that we introduce /3(r) as a local 
r-dependent quantity for reasons which become obvious 
below. To find a reasonable value for /3, which determines 
the local correlation length, we estimate the area where 
the electron is correlated as 



A(r) = d^s e 



-I3\r)s^ _ 



P^r) 



(3) 



Then, we assume that this area is proportional to the area 
of the Wigner circle Trr^ , where the density parameter Ts 
is given through the total electron density as rs(r) = 
l/^TTp{r). Thus, we find the relation 



^(r) =<?VpW 



(4) 



where q is a fitting parameter. 

The SD wavefunction in Eq. ([T]) is recovered when all 
pairs of electrons are far apart from each other. In con- 
trast, when two electrons are brought to the same point, 
the parameter a is chosen to satisfy the cusp condition 



(for the singlet case) of the wavefunction. It can be 
showrJS. that, for the 2D case, a = 1. The function 
controlling the exponential decay is given by 



<j>(r) = 



/3(r) 



/3(r)+0F/2' 
which can be deduced by imposing the conditio n^^i^^ 

d'^s ip{r,s) = 0. 



(5) 



(6) 



By using the wavefunction ([T|) and the definition of the 
correlation factor ip given by ^, we can obtain a formula 
for the correlation energ y ^^i^^ 



Ec = d"^ r d'^s p2,su{r,s) 



(p'^{r, s) — 2ip{r, s) 



, (7) 



where /92,sD(r, s) refers to the SD pair density. To sim- 
plify this expression, we write a Gaussian approximation 
for this function, 



P2,sd(i", s) = P2,SD(l")e 



-sVl'ir) 



(8) 



The use of this Gaussian approximation was proposed 
in the context of the CS scheme by Moscardo and San- 
Fabian,— but it has been used in the field of DFT even 
further backi^ To obtain the function 7(r), that defines 
the width of the Gaussian, we apply the exact sum-rule 



/3SD(r) = T7~T M^^ P2,SD(r,s) 

■p2,SD(r)7'(r), (9) 



iV- 
27r 
N - 1 



from which follows 



7(r) 



'(A^-l)psD(r) 



27rp2,SD(r) 
To simplify this expression, we apply the relation 

P2,SD(r) = ^pioW, 



(10) 



(11) 



as well as Eq. (JH) for the SD density psD(r), and find 
1 



where 



7^(r) 



c = 



c/3^(r), 



2(iV-l)g2 



(12) 



(13) 



Using these results in Eq. ([7]), and performing the inte- 
gration over s, leads to the final result 



i?' 



d^rpsD(r)er"(r). 



(14) 



where we have defined ec(r) as the local correlation en- 
ergy per electron having the lengthy expression 

. (r) = ^|^^Mf$(r) - If + iMM£)^ 



/^$2(j.) 



4/3(r)(2 + c)3/2 yr^ 



^^('•^^^(r)-!]' *^^) 



1 + c 



(15) 



Up to this point, the only inputs for the correlation 
energy are the fitting parameter q (we will come back to 
the choice of this parameter later on), the total number 
of electrons TV, and the electron density p{r). We re- 
mind that the parameter c is defined through q and N 
in Eq. (fT3|) . and (3{r) is given in terms of p(r) in Eq. ([4]). 
This particular depencency on N conflicts with the extent 
requirement of the correlation functional. For example, 
situations where two systems are very far apart from each 
other are expected to be problematic. 

In conclusion, equation (|14p is an explicit density func- 
tional for the correlation energy with a single fitting pa- 
rameter q. This functional is self-interaction free, in the 
sense that it is identically zero for one-electron systems. 
Note that to recover this important property within the 
standard ladder of exchange-correlation functionals, one 
has to resort to highly sophisticated orbital functionals. 



III. APPLICATION AND REFINEMENT OF 
THE APPROXIMATION 

Here we complete and apply the approximation for the 
correlation energy in 2D. In particular, along the appli- 
cations, we shall present an ad-hoc modification which 
better accounts for both the long-range and the kinetic- 
energy correlation. 

As the first step, now we need to choose a value for the 
fitting parameter q, we use Taut's analytic result'^^ for the 
singlet state of a two-clcctron parabolic QD with confin- 
ing strength to = 1. In terms of energy components, the 
correlation energy can be written as Ec = i?tot — EfJ^'^, 
where EXX refers to the exact- exchange result. Applying 
this formula yields Ec « —0.1619 for the TV = 2 singlet 
when oj = 1. To obtain the same value from Eq. (HU, 
we need to set q — 2.258. Of course the choice may be 
refined, if needed. But aiming at providing ideally a pre- 
dictive approximation, the fitting should not be carried 
out for each new system (to obtain every time the cor- 
rect answer) but rather carried out once for ever. This 
is a quite general, and a well known way of defining, or 
refining, new approximations for the central quantities 
of DFT. In the following, we will show that our fitting 
procedure, outlined just above, guarantees a very good 
performance for a large class of systems. 

Tables HI and HIl show results for parabolically confined, 
and for square (tt x n) quantum dots. The results ob- 
tained with our local formula for the correlation energy 



TABLE I: Comparison of the correlation energies (in atomic 
units) for parabolic quantum dots. The reference value E^"^ 
is obtained by subtracting the exact-exchange energy from 
accurate data for the total energy. The last row contains the 
mean percentual error, A, for the parabolic dots (excluding 
the one used in the fitting procedure). 
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0.3805 


0.4453 
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0.3946 
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25.636^ 


26.5528 


0.9168 


0.6837 


0.8504 


1.0000 
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26.1% 


5.9% 


18.4% 



* Fitted result (see text). * Analytic solution by Taut from 
Ref. m. * CI data from Ref. 3. § Variational QMC data 
from Ref. Q. ^ Difltusion QMC data from Ref. g 



TABLE II: Comparison of the correlation energies (in atomic 
units) for square (tt x tt) quantum dots. The reference value 
_Ec°' is obtained by subtracting the exact-exchange energy 
from the quantum Monte Carlo result for the total energy 
(Ref. 135). 



iV E, 



QMC 



-C/t, 



-El 



of 



-£' 



2 3.2731 
6 26.9679 
8 46.7940 



3.4643 0.1908 0.1905 
27.5928 0.6249 0.6578 
47.5962 0.8022 0.9168 



12 103.3378 104.5620 1.2242 1.4494 
16 178.5034 179.9804 1.4770 2.0096 



0.1763 0.2226 

0.5763 0.7624 

0.7836 1.0514 

1.2026 1.6419 

1.6282 2.2534 



A 



19.3% 6.8% 33.6% 



(denoted by £^'°'^^') are compared to reference results 
i?™^, as well as with LDA correlation energies E^^^. We 
computed the EXX and LDA values using the real-space 
code Oct opus. ''^ The EXX result was calculated in the 
Krieger-Li-Iafrate (KLI) approacb24, which is an accu- 
rate approximation in the static case.— For E^'^^ we 
applied the parametrization of Attaccalite et aL— Note 
that wc used a perturbative approach to calculate E}!"^^^ 
from Eq. (fT4|) : The self-consistent EXX density was the 
input for our local functional. We also found that using 
the LDA density as input did not make a considerable 
difference. 

The QDs studied here span a wide range of den- 
sity parameters Tg determined in a parabolic QD as 
Ts = Af~i/Sti;~^/'^ and in our square QD as r^ — ^-k jN . 
This parameter corresponds to the average radius of 
an electron in a QD with an average number density 
TT-o — !/(■'''''?)• Thus, the cases shown in the tables are be- 
tween 0.44 <rs < 9.71. In fact, the upper limit exceeds 
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FIG. 1: Correlation energy per unit particle for the uniform 
2D electron gas in various approximations. 



the threshold of r^ ~ 7.5 for Wigner crystaUization in 
the impurity-containing 2D electron gas^S One should, 
however, bear in mind that in QDs the concept of Wigner 
localization is ambiguous, and no general formula exists 
for the density parameter at the onset of localization. It 
can be also seen that in our examples the ratio of the 
correlation to the total energy varies from less than one 
percent up to around 30%. 

Results with our local formula are roughly of the same 
quality as the LDA, slightly worse for the parabolic dots, 
but slightly better for the square dots. Furthermore, 
a calculation for the homogeneous electron gas (Fig. [T]) 
reveals that this functional agrees with the LDA (which 
is exact for this system) in the limit of vanishing r^ , but 
underestimates the correlation energy otherwise. 

The derived functional not only gives already very rea- 
sonable results, but is also a very good starting point 
for further developments. In fact, we found an alterna- 
tive functional (that we will denote by i?c°mod) obtained 
by modifying the first term between the parenthesis in 
Eq. HSl) by 



[$(r) - if -^ $(r) - 1 



(16) 



To finish the derivation of this new functional we need 
to refit the parameter q, which now reads g™°^ = 3.9274 
(in such a way, we again obtain the exact value of the 
correlation energy for the two-electron quantum dot as 
described above). 

Sensu stricto, Eq. (|16p is an empirical approximation. 
However, our results suggest, Post-factum, that the pro- 
posed modification better accounts for the long-rag o^^i^^ 
and kinetic-energy correlationi22i21 

According to Tables U and [TTl our corrected functional 
agrees very well with the reference results. We find that, 
in all the cases studied, our approximation is vastly su- 
perior to the LDA correlation. Note that our results ex- 
hibit the correct scaling with respect to both confinement 
strength and number of electrons, even if the adjustable 



parameter q has only been fitted to the case A^ = 1 and 
uj = 2. Also for the homogeneous electron gas (Fig. [1]), 
our modified functional yields results that are remarkably 
close to the reference LDA curve, departing significantly 
from the exact curve only for very small r^ (weak corre- 
lation limit). 

Finally, we wish to make a few remarks on the us- 
age of the present correlation functional. First, we point 
out that in practical purposes within, e.g., the Kohn- 
Sham scheme of DFT, the functional should be com- 
bined with an adequate recipe for the exchange energy, 
such as the exact-exchange or the functional suggested 
in Ref. [sl. Second, for many systems — like, e.g., QDs 
in magnetic fields — one requires a spin-polarized version 
of the exchange-correlation functional. This has already 
been taken into account in the LDA functional by At- 
taccalite et alM, but a spin-polarized extension of the 
present functional is still missing. Work to solve these 
two issues is already under way. 



IV. CONCLUSIONS 

We developed a correlation energy functional for the 
two-dimensional electron gas, starting from the CoUe and 
Savetti ansatz for the many-body wavefunction and a 
Gaussian approximation to the pair density. To better 
account for the long-range and kinetic-energy correlation, 
we have then introduced an additional ad-hoc modifica- 
tion. The resulting functional has a very simple form, 
depending parametrically on the total number of elec- 
trons N and locally on the electronic density n(r). It 
only contains a single parameter, q that was adjusted 
to the exact calculation of a two-electron quantum dot. 
Calculations performed for several systems, with a wide 
range of density parameters r^ , show that our functional 
gives results in very good agreement with reference val- 
ues. This agreement is maintained even for very dilute 
electron gases, where the correlation energy amounts to 
30% of the total energy. Furthermore, our functional per- 
forms significantly better than the standard LDA correla- 
tion functional, while maintaining much of its simplicity. 
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